The XMM-Newton wide-field survey in the COSMOS field: V. 
Angular Clustering of the X-ray Point Sources 
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ABSTRACT 

We present the first results of the measurements of angular auto-correlation functions (ACFs) 

of X-ray point sources detected in the XMM-Newton observations of the ~ 2 deg 2 COSMOS field 

(XMM-COSMOS). A significant positive signals have been detected in the 0.5-2 (SFT) band, in 

the angle range of 0.5-24 arcminutes, while the positive signals were at the ~2 and ~ 3<r levels in 

the 2-4.5 (MED) and 4.5-10 (UHD) keV bands respectively. Correctly taking integral constraints 

into account is a major limitation in interpreting our results. With power-law fits to the ACFs 

, „//. 

without the integral constraint term, we find correlation lengths of 9 C = 1.9 ± 0'.'3, 0. 8 ~j, and 

6" ± 2" for the SFT, MED, and UHD bands respectively for 7 = 1.8. The inferred comoving 
correlation lengths, also taking into account the bias by the source merging due to XMM-Newton 
PSF, are r c w 9.8 ± 0.7, 5.8tl\r and 12 ± 2 h^ 1 Mpc at the effective redshifts of z e s « 1.1, 0.9, 
and 0.6 for the SFT, MED, and UHD bands respectively. If we include the integral constraint 
term in the fitting process, assuming that the power-law extends to the scale length of the entire 
XMM-COSMOS field, the correlation lengths become larger by ~ 20%-90%. Comparing the 
inferred rms fluctuations of the spatial distribution of AGNs <T8,agn with those of the underlying 
mass, the bias parameters of the X-ray source clustering at these effective redshifts are in the 
range 6 A gn = 1-5 - 4. 



Subject headings: galaxies: active — (galaxies:) quasars: general 
mology:) large-scale structure of universe — X-rays: galaxies 



cosmology: observations — (cos- 



1. Introduction 

Results from recent X-ray surveys have made 
very significant contributions to understanding 
formation and evolution of supermassive black- 
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holes (SMBHs) at galaxy centers. In particular, 
studies of X-ray luminosity function and its evolu- 
tion have been providing the most reliable current 
estimates of the accretion history to SMBH. One 
of the most important findings in recent years is 
that luminous active galactic nuclei (AGNs) arise 
earlier in the hi story of the universe than lower 
luminosity ones ( Ueda et aD 2003t Hasinger et al 



20051 : iBarger et al.l 120051 lLa Franca et all l2005h . 
This suggests that more massive SMBHs have 
been formed earlier in the universe, and reside qui- 
escently at the centers of giant elliptical galaxies in 
the later epochs, while more numerous, less mas- 
sive SMBHs have been formed and accreted later 
in the history of the universe. 

Clustering properties of AGNs and their evo- 



1 



lution with redshift provide yet additional clues 
to understanding the accretion processes onto the 
SMBHs. These give clues to environments of AGN 
activities. In the framework of the Cold Dark 
Matter (CDM) structure formation scenario, clus- 
tering properties or the bias of AGNs over a suffi- 
ciently large scale 6 A gn = (Sp/p) A GN/(Sp/p) 
may be related to the typical mass of dark mat- 
ter halos in which the y reside, (|Mo fc White! 
19961: ISheth et al] l200ll ). At the same time, 
the mechanisms of triggering the AGN activ- 
ity, which might be closely related to galax> 



interactions and/or m erging (jMenci et al 



galaxy 
.1 \2004 . 



Pi Matteo et al.ll2005l ). yield a clustering of AGNs 



and can therefore be infered from the clustering 
analysis. 

Since strong X-ray emission is a typical fea- 
ture of an AGN activity, X-ray surveys provide 
most efficient means of constructing comprehen- 
sive complete samples of AGNs without contami- 
nation from the light in the stellar population of 
host galaxies. In particular, surveys in the harder 
(E > 2 keV) X-ray band such as available from 
XMM-Newton are very efficient in finding not only 
unobscured AGNs, which are relatively easy to se- 
lect also in the optical bands, but also obscured 
ones, which are difficult to select with optical se- 
lection criteria alone. This is im portant because 
most of the accretion { 



80%; 



Com astri et al 

19951 : iGilli et al.l 120011 : lUeda et alTl2003h occurs 



in AGNs obscured by gas (in X-ray bands) and 
dust (in the optical bands). While one approach 
in investigating the environment of AGNs is to 
measure AG N overdensities around known clusters 
of galaxies (ICappi et al.ll200ll : iD'Elia et all 12004 
ICappelluti et al.l 20051 ). a more common and di- 
rect measure can be obtained by calculating auto- 
correlation functions (ACFs) of well-defined sam- 
ples of X-ray selected AGNs. 

While small number statistics limits the accu- 
racy of the clustering measurements of X-ray se- 
lected AGNs, there are a number of reports on 
the detection of the correlation signals. Samples 
based on the ROSAT All-Sky survey have mainly 
constrained the clustering properties of type 1 lo- 
cal AGNs at z < 0.3. Th e correlation length s re- 
sulti ng from the angul ar (jAkvlas et al. 2000) and 
3D (jMullis et al.ll2004r) analyses of these samples 



are 6-7 hr 1 Mpc 0. Due to the wide redshift 
distribution, it is more difficult to obtain cluster- 
ing signal in deeper X-ray surveys before redshifts 
for a complet e set of X-ray sources a re obt ained. 
Nevertheless, iBasilakos et al.l ()2004 120051 ) mea- 
sured strong angular correlation signals in their 
XMM-JVewton/2dF survey, which covers a total 
area of 2 deg 2 over two fields. They obtained 
a correlation length of ~ 7.5 hr 1 Mpc in phys- 
ical units for both soft and hard X-ray selected 
sources, suggesting a clustering evolution which is 
fixed in the proper coordinate betwee n and 
z > 1. At much fainter X-ray fluxes, IGilli et al 



(2005) analyzed the projected-distance correlation 
function w(r p ) for the X-ray sources with spectro- 
scopic redshifts in the Chandra Deep Fields North 
(CDF-N) and South (CDF-S). They found signif- 
icantly different clustering properties in these two 
fields, suggesting a cosm ic variance effect. Re- 
cently, lYang et al (|2006l) made detailed analysis 
on their 0.4 deg 2 Chandra Large Area Synoptic X- 
ray Survey (CLASXS) supplemented by CDF-N. 
With spectroscopic redshifts for a good portion of 
the sources, they explored the clustering proper- 
ties in different redshift and luminosity bins as well 
as intrinsic absorption bins. They found the evo- 
lution of bias with redshift but they did not find 
significant dependence in the clustering properties 
of X-ray selected AGNs based on either luminosity 
or intrinsic absorption. 



One of the main aims of the COSMOS ( Scoville et al 



2007) project is to trace the evolution of the large- 



scale structure of the universe with an unprece- 
dented accuracy and redshift baseline. The XMM- 
Newton Survey, coveri ng the entire COSM OS 
field (XMM-COSMOS; lHasinger et ail |2007|) . is 
one of the most extensive XMM-Newton Sur- 
vey programs conducted so far. In the first-year 
XMM-Newton observations, about 1400 X-ray 
point sources have bee n detected and cataloged 
(|Cappelluti et all 120071 ) (hereafter C07), which 
are dominated by AGNs at redshifts 0.7 < z < 2 
(|Brusa et alj|2007t iTrump et al.ll2007l) . 

In this paper, we report the first results of our 
investigations on the large scale structure through 
an angular auto-correlation (ACF) function anal- 
ysis of the X-ray point sources detected in XMM- 



1 Throughout this paper we adopt (_ffo,f^M,^A) 
(100 hkms- 1 Mpc" 1 , 0.3, 0.7) 
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COSMOS, as a preview of more detailed studies 
in the near future. Our future studies include the 
derivation of the direct three-dimensional corre- 
lation function using redshift information already 
available for a large portion of the X-ray sources 
and the analysis of the cross-correlation of the X- 
ray sources with galaxies in the multiwavelength 
COSMOS catalog. 

The outline of the paper is as follows. In Sect. 
[51 we explain the selection of our samples of X-ray 
sources to be used in the correlation function anal- 
ysis, which are subsets of those described in C07. 
Details of the calculations, including the ACF es- 
timator, the random sample, and power-law fits 
are presented in Sect. [31 The de-projection to 
the three-dimensional correlation function is pre- 
sented in Sect. [U The results are discussed in 
Sect. [SJ We summarize our conclusions in Sect [HI 

2. Sample Selection 

Our samples consist of the X-ray sources de- 
tected in the first-year XMM-Newton observa- 
tions of the COSMOS field. The source detection, 
construction of the sensitivity maps, and source 
counts are described in C06. The X-ray source 
catalogs in three energy bands, corresponding to 
energy channels of 0.5-2 (SFT), 2-4.5 (MED) and 
4.5-10 (UHD) keV are used. 

For the angular ACF studies, we have applied 
further selection criteria to the C06 sources to 
minimize the effects of possible systematic errors 
in the sensitivity maps. The applied criteria for 
this kind of analysis should be stricter than those 
adopted for the derivation of the log A — log S 
function, because localized systematic errors may 
cause spurious clustering of X-ray sources. In or- 
der to do this, we have scaled up the original sen- 
sitivity map to: 



CRk 



max(a Ci?n miC 06 - b, CR\ 



lm.mm ) > 



(1) 

where CRn m ,C06 is the limiting count rate ( in 
cts s _1 ) in the original C06 sensitivity map and 
C-Riim.acf is the sensitivity map used for the ACF 
analysis. After a number of trials, the scaling co- 
efficients (a, b) have been set to (1.33,1 x 10~ 4 ), 
(1.40, 0.) and (1.44, 4 x 10~ 4 ) for the SFT, MED 
and UHD bands respectively. We have excluded 
the area where CRy lm , ac [ exceeds Ci?ii mjmax (low 
exposure areas close to the field borders). Those 



X-ray sources with Ci?'s below Ci?ii mjac f at the 
source position have been excluded from the ACF 
analysis. After these screenings, the numbers of 
sources for the ACF analysis are 1037, 545, and 
151 for the SFT, MED, and UHD bands respec- 
tively. While the sensitivity in the soft band 
is the best among the three bands, some X-ray 
sources are hard enough that they are detected in 
only MED and/or UHD bands. These are mainly 
highly obscured AGNs. Out of the 545 (151) MED 
(UHD) band sources after the screening process, 
59 (13) have not been detected in the SFT band, 
and only one UHD sources have escaped the detec- 
tion in the MED band (before the screening pro- 
cess). The numbers of the X-ray sources, values of 
Ci?ii m ,min,Ci?iim,max and the total areas used for 
the ACF analysis are summarized in Table [TJ 

3. Angular Correlation Function Calcula- 
tion 

3.1. The ACF calculation 

In calculating the bin ned ACF, we have use d 
the standard estimator by lLandv fc Szalavl (|1993i ): 



Wcst (Qi ) = (DD - 2DR + RR) /RR, (2) 

where DD, DR, and RR are the normalized num- 
bers of pairs in the i-th angular bin for the data- 
data, data-random, and random-random samples 
respectively. Also we use the symbols D and 
R to represent the data and random samples 
respectively. Expressing the actual numbers of 
pairs in these three combinations as n pa i r) DD(#i)) 
n pa ir,DR(0i) and n pa i r) RR,(0j), the normalized pairs 
are expressed by: 

DD = n pah , DD (0 i )/[^D(A r D - 1)] 
DR = n pah , D R(^)/(iV D A R ) 
RR = n p & t B R (0 i )/[N B .(N Ii - 1)] (3) 

where Ad , and Ar are the numbers of sources in 
the data and random samples respectively. The 
number of objects in the random sample has been 
set to 20 times of that in the data sample. This 
makes the variance of the second and third terms 
of Eq. [2] negligible in the error budget of w cst . 

Our XMM-Newton observations have varying 
sensitivity over the field. In order to create a 
random sample, which takes the inhomogeneity of 
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the sensitivity over the field into account, we have 
taken the following steps. 

1. Make a random sample composed of Nn ob- 
jects, where N-r is an integer times JVrj. 

2. For each random object, assign a count rate 
from a source from the data sample. The 
assignments are made in sequence so that 
the CR distributions of the random and the 
data sample objects are exactly the same. 

3. For each random object, assign a random 
position in the field. If the sensitivity-map 
value at this position (Ci?ii m ,acf from Eq. [TJ 
in units of counts s" 1 ) is larger than the as- 
signed CR, find a different position. Repeat 
this until the position is sensitive enough to 
detect a point-source with the assigned CR. 

As a check on this procedure, we have also cal- 
culated ACFs using two other methods of generat- 
ing random samples. The second method is to as- 
sign the count rates to the random sources drawn 
from a logiV - log S relation (e.g. iMoretti et al 



2003), instead of copying the count rates of the 
data sample. Then the source is placed at a 
random position in the field. If the sensitivity 
limit at this position is higher than the assigned 
CR, this random source is rejected. Another, but 
more sophisticated and computationally demand- 
ing method is to generate random sources based 
on the log N - log S relation, down to a flux level 
much lower than the sensitivity limit of our obser- 
vations. These sources are then fed into a simula- 
tor, taking into account XMM-Newton instrumen- 
tal effects, including position-dependent PSF, ex- 
posure maps, and particle background. The entire 
first-year XMM-COSMOS image has been simu- 
lated and the same source detection procedure as 
that applied to the actual data has been applied 
on the simulated data. A random sample R is gen- 
erated from 10 simulated XMM-COSMOS fields. 
Using the above two methods for random sample 
generations did not alter the results significantly. 
In the following analysis, we show the results ob- 
tained by the first method. 

3.2. Error Estimation and Covariance Ma- 
trix 

We have estimated the errors using the variance 
of the calculated ACF by replacing D in Eq. [2] by 



random samples. A random sample with the same 
number of objects and the same set of count rates 
as D has been drawn independently from R in Eq. 
[3J We denote the random sample as a replacement 
of D during the error search by R 1 to distinguish 
from R. For each angular bin, a 1 a error has 
been calculated as a standard deviation from re- 
sulting ACFs calculated from A run =80 different 
R' samples, which is then multiplied by a scaling 
factor [1 + w{8 i )] 1 / 2 (hereafter referred to as "the 
scaled random errors"). This scaling factor cor- 
rects for the difference between the errors in the 
null-hypothesis case, obtained from R's, and those 
in the presence of the correlation signal. This is in 
line with the observation that the error in each bin 
of the ACF calculated using the lLandv fc Szalav 
(|l993l ) estimator is approximately a Poisson fluc- 
tuation of the number of data-data pairs in the 
bin, i.e., a ~ [1 + iy(0)]/ v /n pa ir,DD- The standard 
deviation of the null- hypothesis ACF, obtained by 
replacing D by an R', is cr ran ~ ■ 
The scaling factor can be obtained by using the 
relation, n pa ir,DD ~ [1 + w(9)]n pa , h . iR , R , . 

In correlation functions, the errors in different 
angular bins are not independent from one another 
and correlations among the errors have to be taken 
into account when we make function fits. Thus we 
have also estimated full covariance matrix in order 
to represent the correlations among errors by 

M coy<i j = J2k( w R,k(0i) - (ws,(6i))) 

x(l + U ; cst (^)) 1 /2(l + u , cst (^))i/2 (4) 

where u>n,,fc(#j) is the ACF value for the fc-th ran- 
dom run (k runs through 1 to N Tun — 80), (wn(9i)) 
is their mean value, evaluated at the center of the 
?-th angular bin 0j and w est (9i) is from Eq. [21 The 
square roots of the diagonal elements of Eq. 2] are 
the scaled random errors discussed above. The co- 
variance matrix calculated in Eq. 0]is used later in 
Sect. 13.41 Strictly speaking, Eq. [4] only takes into 
account the correlations of errors for the random 
cases, but not the correlation of errors due to clus- 
tering. One way of explaining this is that remov- 
ing/adding one source (by a Poisson chance) af- 
fects multiple angular bins and this is represented 
by non-diagonal elements of Eq. [4] On the other 
hand, correlation of errors due to large scale struc- 
ture or the cosmic variance is not represented by 
this. If we observe another part of the sky, we 
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sample different sets of large scale structures such 
as filaments and voids. Since the existence or non- 
existence of one such structure affects multiple an- 
gular bins, there should be additional contribution 
to the non-diagonal elements of the correlations 
among errors in different angular bins. Eq. [4J 
based on many random samples, thus includes the 
former type of the correlation of errors but not 
the latter. One way to include also the latter ef- 
fect is to use the Jack knife re-samp l ing te chnique, 



as was done by e.g. IZehavi et al.1 (|2004f ) . How 



ever, the Jackknife re-sampling requires to divide 
the sample into many statistically independent re- 
gions, which is not practically possible in our case. 

3.3. The binned ACF results 

The ACFs have been calculated for the three 
bands in logarithmically equally spaced bins with 
Alog6> = 1/6. The results are shown in Fig [IJ 
where the upper panels, composed of two layers of 
logarithmic plots with positive and negative parts 
(log|ui(#)| > —2.8 respectively), are attached to- 
gether. The lower panels show fit residuals for the 
best- fit functions described in Sect. 13.41 Changing 
the bin size did not change the clustering ampli- 
tude significantly. 

The ACFs are presented with the scaled ran- 
dom errors. Positive signals have been detected 
down to ~ 0.5' in the 0.5-2 keV band and 6 ~ 1' 
in the other bands. At the smallest scales, cor- 
relation signal is negative, probably due to con- 
fusion effects, where two sources separated by a 
distance comparable to or closer than the point 
spread function (PSF) cannot be detected sepa- 
rately in the source detection procedure and may 
well be classified as one extended source. In our 
current sample, the sources that have been classi- 
fied as extended have been removed from the sam- 
ple. The effect of this is discussed in detail in Sect. 
13.51 In the 0.5-2 and 2-4.5 keV bands, positive sig- 
nals extend out to ~ 20'. Negative signals are seen 
at the largest angular scales, probably due to the 
integral constraint as discussed below. 

3.4. Power-law Fits 

In order to make a simple characterization of 
our ACF results, we have fitted the ACF with a 
power-law model of the form: 



where 7 is the slope index of the corresponding 
three-dimensional correlation function. We use 
the normalization A as a fitting parameter rather 
than the angular correlation length 9 C = A 1 '" -1 ', 
since this gives much better convergence of the fit. 

The fits are made by minimizing \c- The sub- 
script c denotes that the correlations between er- 
rors have been taken into account through the in- 
verse of the covariance matrix: 



Xc 



A T m: 



1 A 



(6) 



where A is a vector composed of w cs t{0i) — 
ifmdi(^i) + C, Mqov is the covariance matrix cal- 
culated in Eq. 2] with N lna = 80, and C is a 
constant to compensate for the integral constraint 
as discussed below. 

Due to the finite area and the construction 
of West in Eq. [21 the estimated angular correla- 
tion function satisfies the integral constra int (e.g 
Basilakos etaHl2005l: iRoche fe Ealeslll999t ): 



w cst d ft = 0. 



(7) 



This constraint usually results in w cs t underes- 
timating the true underlying angular correlation 
function by the constant C. Under an assumption 
that the true underlying w{0) is a power-law and 
is extended to the scale of the survey area, one can 
include C in the fitting process, where C can be 
uniquely determine d by 9 C and 7 by imp osing the 



integral constraint ijRoche fe Eale ; 



by imp os 
3sl ll999h . 



(8) 



Wmdl 



(e) = a e 1 - 



(•5) 



where the sums are over angular bins and RR(0i) 
is the number of random-random pairs in the i-th 
angular bin. The above assumption is not nec- 
essarily true. If residual systematic errors in the 
sensitivity maps are the main cause of the negative 
values at large angular separations, the determina- 
tion of C shown above is not valid. However, this 
gives an approximate estimate of the degree of the 
underestimation by the integral constraint. This 
sets an limitation to the our angular ACF analy- 
sis, where the estimated C values are not negligible 
compared with the amplitude of the ACF signal. 
We have made fits with or without including the 
integral constraint. 

Because of the limited signal-to-noise ratio, we 
were not able to constrain A and 7 simultaneously. 
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Fig. 1. — The binned estimated angular correlation functions w cs t(0) are plotted for the X-ray sources 
detected in the first-year XMM-Newton data in three standard energy bands as labeled. The vertical scale is 
logarithmic, where the positive and negative parts (log |iu e st(#)| > —2.8 respectively) are attached together. 
The ler errors are the diagonal components of Eq. [|J i.e., the scaled random errors. The blue-solid and 
red-dotted lines (colors in the electronic version only) show the best-fit power-law models for 7 — 1 = 0.8 
without and with an integral constraint respectively. The models are plotted in the range where the fits are 
made. Fit residuals in terms of a has been also plotted in the lower panels for the two models (slightly offset 
for visibility) in the same line styles (colors) as the models. 
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Thus we have calculated the best-fit values and la- 
confidence errors for the amplitude for two fixed 
values of 7 — 1 = 0.8 and 0.5. The former value 
is for the canonical value for lo cal galaxies (e.g. 
Peebled[l980t IZehavi et alj|2004 . the latter is ap- 
proximately the slope fou nd for X-ray sou rces in 
the Chandra Deep Fields (|Gilli et al-lboOSl ). The 
angular fit results are summarized in Table [2l In 
this table, fits with different bands and parameters 
are identified with a Fit ID. The angle range for 
the fits are 9 min < 6 < 6* max and the boundaries 
are also shown in Tabled For fit ID's S1-S4, MI- 
AM, and U1-U4, we have set d min = 0.5' (SFT) or 
0.7' (MED,UHD), which is the minimum at which 
ACF is still positive, and below which the ACF 
goes negative due to the XMM-Newton PSF. Like- 
wise, we set # max = 24', which is about the max- 
imum scale where the ACF is still positive. The 
best-fit models for 7 — 1 = 0.8 are overplotted in 
Fig. [T]in the bin ranges included in the fits. 

As another choice, we have set 9 m i n and 6 max 
in such a way that the range approximately corre- 
sponds to the projected comoving distance range 
of 1-16 h- 1 Mpc (Fit ID S5,S6,M5,M6,U5, and 
U6) at the effective median redshift of the sam- 
ple (z e ff discussed below in Sect. 0J The rationale 
for the maximum scale is that, in our subsequent 
analysis, the correlation functions are converted 
to the root mean square (rms) density fluctuation 
with a 8/i _1 Mpc-radius sphere (therefore the rel- 
evant maximum separation is 16ft. -1 ) in discussing 
bias parameters. The rationale for the minimum 
scale is to minimize the effects of non-linearity in 
discussing the bias parameters and typical halo 
masses. 

3.5. Effects of Source Merging due to PSF 

The amplification bias, due to which the esti- 
mated ACF from sources detected in a smoothed 
image (e.g. by a finite PSF) is amplified with 
respect to the true underl ying ACF, has been 
first n oted and discussed by lVikhlinin fc Formanl 
(|1995l ) in the context of the clustering of X-ray 
sources. This is caused by merging of multiple 
sources which are separated by distances compa- 
rable to or closer than the PSF. The effect of this 
bias depends on the true underlying angular cor- 
relation function and the number density of the 
sources. In principle, full simulations involving 
PSF smoothing and the source detection process 



are required to estima te the amount of this bias. 



Basilakos et al.l (|2005f ) took a simplified approach 



in estimating this effect on their ACF from their 
XMM-JVewton/2dF survey. In order to estimate 
the size of the effect, they used particles sampled 
from a cosmological simulation. They simulated 
the XMM-Newton sources by merging all the par- 
ticle pairs closer than 6". They then compared 
the angular ACFs from the particles themselves 
and the simulated XMM-Newton sources. As a 
result, they estimated that the measured angular 
correlation length is overestimated by 3-4% due to 
the amplification bias. 

In our case with XMM- COSMOS, we have ex- 
plicitly excluded sources that are classified as ex- 
tended by the source detection procedure (C06). 
This causes most of the source pairs closer than 
~ 20" to disappear from the sample, since these 
pairs are classified as single extended sources. Be- 
cause the exclusion of these sources can suppress 
the estimated angular correlation function, we use 
the term "PSF merging bias" rather than the "am- 
plification bias". Pairs of sources that are closer 
than ~ 4 arcseconds are, however, detected as a 
single poi nt source. We have ap plied a similar ap- 
proach to lBasilakos et al.l (|2005l ) in estimating the 
effects of the PSF merging bias. We have sam- 
pled particles from the COSMOS-Mock catalog 
extracted from the Millennium simulation (Kit- 
bichler, M., priv. comm.) over ~ 2 deg 2 of 
the sky. Redshift, cosmological intrinsic redshift, 
and magnitudes in various photometric bands are 
provided for each mock galaxy in the catalog. 
We use the mock catalog to estimate the effect 
of the PSF merging bias in the angular correla- 
tion function. Thus the selected objects from the 
mock catalog for our simulation do not have to 
physically represent to the actual X-ray selected 
AGNs. For our present purpose, we have chosen 
the mock galaxies in a redshift interval (roughly in 
the range 0.4 < z < 0.8) and a magnitude range 
in such a way that the amplitude of the result- 
ing angular ACF and the source number densities 
roughly match those of the X-ray samples. We 
have then created a simulated XMM-COSMOS 
catalog as follows: 1) source pairs with separa- 
tions smaller than 4 arcseconds are merged into 
single sources and 2) pairs that are between 4 to 
20 arcseconds from each other are eliminated. We 
repeated this experiment 19 times and compared 



7 



the mean angular ACFs from the original particles 
and that from simulated XMM-COSMOS by mak- 
ing power-law fits to the mean ACFs. As a result 
we found that the ACF amplitudes measured using 
the sources in our source detection procedure on 
the XMM-COSMOS data are underestimated by 
15% and 8% for the SFT and MED bands respec- 
tively. Corrections for this effect have not been 
applied for the values in Table [21 but are consid- 
ered in further discussions. The effects is negligible 
in the UHD band, due to the relatively low num- 
ber density of the sources detected in this band, 
which made the average distance among neighbor- 
ing sources much larger than the PSF. 

4. Implication for 3-D Correlation Func- 
tion and Bias 

4.1. De-Projection to Real Space Correla- 
tion Function 

The 2-D ACF is a projection of the real-space 3- 
D ACF of the sources £(r) along the line of sight. 
In the following discussions and thereafter, r is 
in comoving coordinates. The relation between 
the 2-D (angular) ACF and the 3- D ACF is ex- 
pressed by the Limber's equation ( Peebles! 1980l 
e.g.,). Under the usual assumption that the scale 
length of the clustering is much smaller than the 
distance to the object, this reduces to: 



i{^[d A {z)9f + P (1 + z)) {^Y 1 dl dz, (9) 



where dx{z) is the angular distance, N is the to- 
tal number of sources and dN/dz is the redshift 
distribution (per z) of the sources. The redshift 
evolution of the 3-D correlation function is cus- 
tomarily expressed by 



S = H. 



£(r,z) = (r/r Cl0 )-i(l + z) 



-3-e+7 



(10) 



where e = — 3 and e = 7 — 3 correspond to the case 
where the correlation length is constant in phys- 
ical and comoving coordinates respectively. In 
these notations, the zero-redshift 3-D correlation 
length r C fi can be related to the angular correla- 
tion length 9 C by: 



J(^) 2 ( £ ^ £l )"^ 7 (l + -)- 3 -^ 

(11) 



„ _ r[( 7 -i)/2]r(i/2) 
n 't ~ r( 7 /2) 



where t(z) is the look back time. Note that all de- 
pendence on cosmological parameters are included 
in dx{z) and t{z). We also define the comoving 
correlation length: 

r c (z eS ) = r Cl0 (l + z cff ) ( - 3 - £+ ^ (12) 

at the effective redshift z e g, which is the median 
redshift of the contribution to the angular corre- 
lation (the integrand of the second of Eq. [TTj) . 
0.8 




Fig. 2. — The model redshift distributions of the 
COSMOS-XMM sources in the 0.5-2 keV (solid 
lines), 2-4.5 keV (dotted line), and 4.5-10 keV 
(dashed line) respectively. The thicker and thin- 
ner solid lines correspond to t he 0.5-2 keV band 
redshift di stributions based on |Ueda et al.1 (|2003f ) 
model and lHasinger et al] ( 20051 ) soft X-rav lumi- 



nosity function respectively. 

An essential ingredient of the de-projection pro- 
cess is the redshift distribution of the sources. At 
this stage, we do not yet have individual redshifts 
of a comprehensive set of the XMM-Newton com- 
plete sample. Thus we use expected distributions 
from the X-ray luminosity functions and AGN 
po pulation synthesis models. We use the model 



po pulation syntnesis models, we use tne model 
bv lUeda et al. ( 2003f) (Luminosity-dependent den- 
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r> ' 



(N 2 /s)er 



sity evolution or LDDE) for all bands. In cal- 
culating the redshift distribution, we have used 
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the sensitivity map in units of CR and the actual 
XMM-N ewton response f unctio n in each band. We 



also use iHasinger et al. (|2005h type 1 AGN soft 
X-ray luminosity function (SXLF) for the 0.5-2 
keV for comparison. The rcdshift distributions 
of the X-ray sources predicte d by these models 
are plotted in Fig B oth llleda et alJ {2003) 



and IHasinger et al.l (|2005h used samples with a 



very high identification completeness with red- 
shifts measurements (> 90%), at least down to the 
flux limits sampled by the XMM-COSMOS sur- 
vey. Thus the effect on the expected redshift dis- 
tribution due to the identification incompleteness 
is negligible. . In calculating the three dimen- 
sional correlation functions, we use the fits with 
and without integral constraints. The angular cor- 
relation amplitude A has been multiplied by a cor- 
rection factor due to PSF merging as discussed in 
Sect. 13.51 Also we use the fits with 7 = 1.8. 
The calculated r c .o and the comoving correlation 
length at the effective median redshift r c (z e g) are 
listed in Table [3] for selected results. The errors 
on r Ci o and r c (z e g) have been calculated for fixed 
7 and e. 

4.2. Bias and Comparison with Other 
Works 

In order to estimate the bias parameter of the 
X-ray sources with respect to the underlying mass 
distribution, we calculate the rms fluctuation of 
the distribution of the X-ray sources in the sphere 
with a comoving radius of = 8/t -1 for the 

power- law model (e.g. §59 of |Peeblej[l980l) . 



°8,agn =/K(|ri-r 2 |)dW 2 /V 2 

= (r max /r c )-^ J 2 (13) 
J 2 =72/[(3- 7 )(4- 7 )(6-7)n (14) 

As discussed above, we have used the results 
from Fit ID's S5, S6, M5, M6, U5,and U6, where 
the fits are made to the angle range correspond- 
ing to » 1-16 hr 1 Mpc at z c g. The correspond- 
ing quantity of the underlying mass distribution 
at z = 0, a s is one of the commonly use d pa- 
rameters in cosmology (jSpergel et al.1 120031 ). In 
order to compare our results with other simi- 
lar works on a common ground, we calculated 
08, agn from power-law representations from lit- 
erature and plotted them versus the effective red- 



shift of each sample^. 

For this comparison, we have used the best fit 
correlation lengths and slopes (r c ,7) from litera- 
ture to estimate <7s,agn values and their la errors. 
Since each reference has a different method of pre- 
senting results, we take the following strategy in 
calculating <7s,agn and its la error. 



(a) 



If the referenced article gives confidence 
contours in the (r c ,7) space, we calculate 
08, agn values for the nominal case as well 
as at each point in the L = L lri \ n + 1 contour, 
where L (with the best-fit value Z/ m in) is ci- 
ther x 2 or a statistical estimator that varies 



as y . e.g. Cash C-estimator for lYang et al 



(|2006l ). The error range on ct"8,agn i s de- 
termined by the minimum and maximum 
values calculated from the points along the 
contour. 



(b) If no confidence contour in the (r c ,7) space 
is given and there is a fit result with a fixed 
7, we use this fit to calculate <js,agn- The 
la error in r c is propagated from that of r c . 

(c) If the article gives only best-fit (r c , 7) values, 
and la errors on both parameters, the error 
of (78, agn has been propagated from those 
of r c and 7, neglecting possible correlation of 
errors between these two parameters. In this 
case, we may well have over/under estimated 
the errors on cts,agn- 

If the (r c , 7) values are given in multiple evo- 
lution models, we use the one where the correla- 
tion length is fixed in the comoving coordinates 
(i.e. e = —3 — 7). At about the effective me- 
dian redshift of the sample, however, the correla- 
tion lengths calculated assuming different values 
of e do not differ significantly. In the case of our 
present work, we see this by comparing the r c (z e s) 
values for e = —1.2 and e = — 3 cases in Table [3] 
Also the value of cts.agn is insensitive to the as- 
sumed value of 7. The change of <7s,agn is less 
than 0.1 between the assumed 7 of 1.8 and 1.5 
for our results in all the three bands. The results 



2 Some authors give the median redshift of the number dis- 
tribution of the X-ray sources, while we and some others 
give median redshift of the contribution to the clustering 
signal. We denote the former by z and the latter by z e ff. 
We do not make a distinction between these in Fig. \3\ 
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from literature we use for this comparison and the 
details of the conversion to o^agn are described 
below, roughly in order of redshift. 

Grazian et al.l ( 20041 ) calculated the correlation 



function of 392 optically-selected QSOs from the 
Asiago-ESO/RASS QSO Survey (AERQS) with 
z cS = 0.062. They found the nominal values of 
r c = 8.6 ± 2.0 hr 1 with 7 = 1.56. Also in the low- 
redshift end, lAkvlas et al. ( 2000h calculated the 
correlation function of X-ray selected AGNs from 
the ROSAT All-Sky Survey (RASS) with a me- 
dian redshift of z — 0.15. Their correlation length 
for 7 = 1.8 of r c = 6.5 ± 1.0 h' 1 Mpc for the 
Einstein de-Sitter Universe is increas ed by 5% to 



conve rt it to our adopted cosmology. iMullis et al 
(<2004l) found r c = 7.4±i;| h' 1 Mpc for 7 = 1.8 in 



their ROSAT North Ecliptic Pole Survey (NEPS) 
AGNs with median redshift for the contribution to 
the clusterin g signal of z c g = 0.22. iBasilakos et al.l 
(|2004l2005l) in their 2 deg 2 XMM-Newton survey, 
with shallower flux limits than XMM- COSMOS, 



found 



= 16.4 ± 1.3 h- 



at z 



— 1.2 and r c = 



19±3 h- 1 at z = 0.75 for th e 0.5-2 and 2-8 keV re - 
spectively. A recent work bv rPuccetti et al. ( 2006h 
on the central ~ 0.6 deg 2 region of the ELAIS-S1 
field , covered by four mosaiced XMM-Newton ex- 
posures with ~ 50 — 60 ks each, also measured 
angular ACFs of X-ray point sources. For fixed 
7 = 1.8, they found r c = 12.8 ± 4.2 h' 1 at z = 1.0 
and r c = 17.9 ± 4.8 h~ l at z = 0.85 for the 0.5-2 
and 2-10 keV bands respectively. 

The correlation functions on the deepest X-ray 
surveys on the Chandra Deep Fields- South (CDF- 
S; z = 0.84) and North (CDF-N; z = 0.96) by 
Gilli et all (|2005h gave, for fixed 7 = 1.4, r c = 
10.4 ± 0.8 h' 1 and 5.1±g;|/i -1 Mpc respectively. 
We use the results from their AGN samples. For 
all of the above samples, the errors on the o^agn 
have been calculated using method (b). 

An extensi ve redshift-space co rrelation function 
was made by lYang et all (120061) . who made use 



of the data from a combination of the CLASXS 
and CDF-N surveys, with a significant portion of 
the X-ray sources having measured spectroscopic 
redshifts. We have used their (so,7) confidence 
contours, where sq is the redshift-space comoving 
correlation length, in the four redshift bins with 
median redshifts of z = 0.45, 0.92, 1.26, and 2.07 
to estimate o^agn using method (a). In this con- 
version, we have corrected for the redshift distor- 



tion by dividing the redshift-space cr s . a on valu e 
by \/r3 (jMarinoni et al.l 120051: lYang et al.ll2006l ). 
Extensive clustering studies of QSOs from the 2dF 
QSO redshift survey (2QZ) have been made using 
both the projected-distance co rrelation function 
approach (|Porciani et al 1 l2004h and the redshift- 
space t hree-dimensional co rrelation function ap- 
proach (|Croom et al. l2005h . We have converted 



the nominal ro-7 values and confidence contours 
in three redshift bins at z e ff=l-06, 1.51 and 1.89 
by IPorciani et all (2004) to as. AON using method 
(a). In converting the lCroom et al.l (|2005l )'s (so, 7) 
results in 10 redshift bins ranging from z = 0.5 to 
2.5, we have used method c) and the redshift dis- 
tortion correction has been made in the s ame way 
as we have done to the Yang et al. ( 20061 ) results. 



Figure |3fa) shows the o^agn values as a func- 
tion of the look back time t{z) for our default 
cosmology from our analysis results both with- 
out and with integral constraints. We also over- 
plot o^agn values calculated from the results 
found in literature as detailed above. In or- 
der to compare them with those of the under- 
lying mass distribution, we have also plotted 
the erg D(z) values from the linear theory (e.g. 
ICarroll et all Il992t Irlamiltonl l200lh normalized 
to 0.75 at z=0 (jSpergel et al.ll2003f )Pl. This curve 
has been shown to accurately represent the dis- 
tribution of dark matter p articles in the ACDM 
Hubble Volume Simulation (jMarinoni et al.ll2005h . 
Figure [3)Jb) shows the inferred bias parameters 
&agn = (T8,AGN(z)/[a 8 D(z)]. The values of 
C8,agn and &agn from this work are shown in 
Table H 

5. Discussion and Prospects 

In this work, we used all the point sources 
above the scaled sensitivity threshold without fur- 
ther classification of the sources. We analyzed 
our results assuming that all the X-ray sources 
are AGNs. This, in practice, is a good approx- 
imation. Our preliminary identifications of the 
sources indicate that out of the 1037, 545, and 151 
sources selected for the ACF analysis for the SFT, 
MED and UHD bands, 20, 5, and 1 are appar- 



3 We use the latest value of ug as of writing this paper 
obtained from 

http: //lambda. gsfc .nasa.gov/product/map/current/parameters . cf m 

for the A-CDM model derived from all datasets. 
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Fig. 3. — (a) The o^agn of the X-ray sources/AGNs inferred by the power-law fits to the correlation 
functions from this work and literature are plotted against the look back time corresponding to the effective 
median redshift of the samples. All error bars indicate la errors. The results from this work for the fits 
without integral constraints are plotted with large solid symbols as labeled (with colors in the electronic 
version). Those with integral constraints are also plotted with smaller symbols and at positions slightly 
shifted rightward for visibility. The o^agn values calculated from various results in literature (see text) 
are also shown as labeled (abbreviated as the first author followed by the last two digits of the publication 
year). The dotted line shows <t$ D(z) for the mass in the linear theory normalized to 0.9 at z=0. (b) The 
bias parameter 6agn = &8,agn{z)/[<T8 D(z)\ are plotted as a function of the effective redshift. The meaning 
of the symbols are the same as panel (a). We also show redshift ticks at the bottom part of the figure. 
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ent Galactic stars respectively. Removing these 
sources from the analysis changed the results very 
little. Also our results are not likely to be heavily 
affected by the contamination of clusters/groups, 
since these sources are extended by > 20" (e.g. 
Finoguenov et al.ll2007l ) and are likely to be clas- 
sified as extended by the source detection proce- 
dure, hence removed from our sample. 

As seen in Fig. [3J with an exception of the MED 
band, our analysis without integral constraints 
gives somewhat larger cts,agn values than those 
obtained from results using 2dF opti cally selected 
QSOs by Porciani et al. ( 20041 ) and Croom et al 



(|2005l ). but in general agreement wit h the values 
from Chandra CL ASXS+CD F-N by [Yang_etal' 



(<2006l) . CDF-S bv iGilli et all (I2005D and XMM- 
Newton results from iPuccetti et al.l ( 20061 ). Most 
likely due to the c osmic variance over a small FOV, 
Gilli et alJ (|2005l) 's result on CDF-N gave a signif- 
icantly smaller correlation amplitude than their 
own CDF-S values as well as our results. The an- 
gular ACFs from a s hallower XM M-Newton survey 
bv iBasilakos et al.l (|2004l l2005h gave significantly 
larger <78,agn values than other works in both 0.5- 
2 keV and 2-8 keV bands. The reason for their 
distinctively large value is unclear. 

One of the interesting questions in investigat- 
ing clustering properties of X-ray selected AGNs 
is to investigate whether there is any difference 
in the environments of obscured and unobscured 
A GNs. Apply i ng th e population synthesis model 
of lUeda et al.l (|2003l ) to our sensitivity maps, only 
~ 17% of the sources detected in the SFT band 
at z ~ 1 are obscured AGNs with N H > 10 22 
cm~ 2 . The fraction increases to ~ 40% in the 
MED and UHD bands. A comparison of bias 
parameters between SFT band and MED band, 
which have similar z e s values, seems to show a 
lager bias parameter for the SFT sample. How- 
ever, with the combination of statistical uncer- 
tainties and uncertainties in modeling the inte- 
gral constraint in the MED band, we can only 
conclude that the bias of the obscured AGNs is 
not stronger than that of unobscured AG Ns. In 



other wo rks. IGilli et al.1 (120051); lYang et all (|2006h 
as well as IBasilakos et al.l (I2004LI2005I ) did not find 



any statistical difference between the clustering 
properties of these two. Further studies involv- 
ing the second-year XMM-COSMOS data, which 
in effect doubles the XMM-Newton exposure over 



the COSMOS field, and redshift information of in- 
dividual objects will probe into this problem fur- 
ther. Also with the accepted C-COSMOS program 
totaling 1.8 Ms of Chandra exposure, we will be 
able to probe the correlation functions to a much 
smaller scale, enabling us to investigate the imme- 
diate neighbor environments of these AGNs. 

Our measured bias parameters based on the 
rms fluctuations in the 8/i -1 Mpc radius sphere are 
in the range 6agn = 1-5 — 4. The clustering prop- 
erties of dark matter halos (DMH) depend on their 
mass (|Mo k White! Il996l ISheth et all l200ll ) and 
we can estimate the typical mass of the DMHs in 
which the population of AGNs represented by our 
sample reside, under the assumption that the typ- 
ical mass halo is the main caus e of the AGN bias- 
ing. Following the appr oach of Yang et al. ( 2006h 



and Croom et al. ( 2005) who utilized the model by 



Sheth et al.l (l200ll ). we roughly estimate that the 
typical mass of DMH is ~ 10 13 - 10 14 M for our 
SFT and UHD samples (see Table H]). These are 
an order of magni t ude la rger than those estimate d 
bv lPorciani et al.1 (|2004l ) and lCroom et alJ (|2005l ). 
probably reflecting the large bias parameters from 
our results. 

One of the largest uncertainties in our analy- 
sis lies in the treatment of the integral constraint, 
because its effect is not negligible in our case com- 
pared with the ACF amplitudes in the range of 
our interest. Fig. [3] shows that our results based 
on the fits with integral constraint, under an as- 
sumption that the fitted power-law behavior of the 
underlying w(0) extends to the scale of the entire 
FOV, give a somewhat larger correlation ampli- 
tudes. This assumption may not be true. Also 
the apparent negative w es t(0) values at 6 > 30 in 
Fig. [1] may well be caused by remaining system- 
atic errors. Thus the interpretation of the angular 
correlation functions, where the signals are diluted 
by the projection along the redshift space, has a 
major limitation in correctly taking the integral 
constraint into account. 

The situation will improve when redshift in- 
formation on individual X-ray sources becomes 
available for a major and comprehensive set of 
the X-ray sources. This will enable us to cal- 
culate three-dimensional correlation functions or 
projected-distance correlation functions in a num- 
ber of redshift bins. With the line-of-sight di- 
lution effect suppressed, we will be able to ob- 
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tain a larger amplitude in the correlation signal, 
making the analysis much less subject to the un- 
certainties in the integral constraint. With opti- 
cal followup programs underway on the COSMOS 
field through Magellan and zCOSM OS projects 
(jTrump et al.ll2007t iLiilv et alj|2006j) . we are ob- 
taining spectroscopic redshifts from a major frac- 
tion of the X-ray sources. At the time of writing 
this paper, we have been able to define a sample 
of 378 XMM-COSMOS detected AGNs with mea- 
sured spectroscopic redshifts (~ 30% of the X-ray 
point sources), with a median redshift of z ~ 1. 
Our preliminary analysis of these sources based on 
the projected-distance correlation function w(r p ) 
gives a comoving correlation length of r c sa 8/i _1 
Mpc and 7 = 1.6, which is fully consistent with 
our results without the integral constraints. The 
results of a full analysis utilizing the redshift infor- 
mation will be presented in a future paper (Gilli 
et al. in preparation). 

Extensive multi-wavelength coverage and the 
availability of a galaxy catalog also enables us to 
investigate the cross-correlation function between 
X-ray selected AGNs and galaxies. By cross- 
correlating the X-ray selected AGNs with three 
orders of magnitude larger number of galaxies, we 
will be able to investigate the environments of the 
AGN activity in various redshifts with much bet- 
ter statistics. 

6. Conclusions 

We have presented the first results on the an- 
gular correlation functions (ACFs) of the X-ray 
selected AGNs from the XMM-COSMOS survey 
and reached the following conclusions. 

1. A significant positive angular clustering sig- 
nals has been detected in the 0.5-2 (SFT) 
bands in the angle range of 0.'5-20', while 
in the 2-4.5 (MED) and 4.5-10 keV (UHD) 
bands, the positive signals are 2 and 3<7 re- 
spectively. The robustness of the estimated 
correlation functions has been verified us- 
ing different methods of generating random 
samples. 

2. Power-law fits to the angular correlation 
function have been made, taking into ac- 
count the correlation of errors. Correctly 
taking the integral constraint into account is 



a major limitation on interpreting the angu- 
lar correlation function. For fits with fixed 
7 — 1 = 0.8 and without (with) the inte- 
gral constraint term, we found correlation 

lengths of 6 C = l'.'9±0'/3, 0'/8 +0 'v 5 and 6"±2" 
b ' -0.4 

(3'.'1 ± 0'.'5, 2'.'2 ± l'.'O, and 14" ± 5") for the 
SFT, MED, and UHD bands respectively. 

3. Due to XMM-Newton PSF, most of the 
source pairs closer than ~ 20" are classi- 
fied as single extended sources, and there- 
fore excluded from the sample. This causes a 
bias in angular correlation function measure- 
ments. We have estimated this effect (the 
PSF merging bias) by simulations and found 
that the estimated ACF underestimates the 
amplitude of the true underlying ACF by 
- 15% and - 8% for the SFT and MED 
bands respectively 

4. Using Limber's equation and the expected 
redshift distributions of the sources, we have 
found comoving correlation lengths of r c 
9.8 ± 0.7, 5.8±^;f, and 12 ± 2 h' 1 Mpc for 
7 = 1.8 at the effective redshifts of z c s w 1.1, 
0.9, and 0.55 for the SFT, MED, and UHD 
bands respectively for the fits without inte- 
gral constraints, while 20%-90% larger cor- 
relation lengths have been obtained for the 
fits with integral constraints. 

5. Using the fits in the angles corresponding to 
a projected distance range of l-16/i _1 Mpc 
at the effective median redshift of the sam- 
ple, we have calculated the rms fluctuations 
of the X-ray source distributions. Compar- 
ing them with that of the mass distribution 
from the linear theory, we find that the bias 
parameters of the X-ray sources are in the 
range &agn = 1.5 — 4 at 0.5 < z < 1.2. 

6. If the bias mainly reflects the typical mass 
of dark matter halos in which these X-ray 
AGNs reside, their typical masses are 10 13 — 
10 14 M Q . 

7. Further investigations utilizing redshifts of 
individual X-ray sources and/or involving 
cross-correlation function with galaxies tak- 
ing advantage of the wealth of multiwave- 
lcngth data are being conducted. The ap- 
proved Chandra observations (C-COSMOS) 
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on this field will enable us to probe into the 
clustering in much smaller scales and there- 
fore into immediately neighboring environ- 
ments of AGNs. 
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Table 1 

x-ray sources and sensitivity limits 



Band 


Number 


^-^lim.min C^lim,max 


S x> u m range 




Area 


(keV) 




(ctss" 1 ) 


(erg s -1 cm -2 ) 


a 


(deg- 2 ) 


SFT 


1037 


7.0 10" 4 -2.2 10~ 2 


6.7 10~ 16 -2.1 10 


-14 


1.43 


MED 


545 


7.0 10~ 4 -2.5 10~ 2 


4.6 10~ 15 -1.6 10 


-13 


1.56 


UHD 


151 


9.0 10" 4 -2.2 lO" 2 


8.7 10~ 15 -1.8 10 


-13 


1.25 



a Flux range corresponding to the CR limits. The conversions have been 
made, following C06, to fluxes in 0.5-2.0, 2.0-10, 5.0-10 keV assuming power- 
law spectra of photon indices V = 2.0, 1.7, and 1.7 for the SFT, MED, and 
UHD bands respectively. 



Table 2 

Results of the Power-law fits 



Fit ID 


Band 




7-1 


9 c a 


C 


$min 


Pmax 




(keV) 


(arcmin 1 /^ 7 - 1 )) 




(") 




(') 


(') 


SI 


SFT 


0.063 ± .008 


0.8 


1.9 ±0.3 





0.5 


24 


S2 


SFT 


0.093 ± .012 


0.8 


3.1 ±0.5 


8 lO" 3 


0.5 


24 


S3 


SFT 


0.034 ± .004 


0.5 


0.07 ± .02 





0.5 


24 


S4 


SFT 


0.078 ± .010 


0.5 


0.37 ± .09 


1.5 10~ 2 


0.5 


24 


S5 


SFT 


0.059 ± .009 


0.8 


1.7 ±0.3 





1.6 


24 


so 


SFT 


0.089 ± .015 


0.8 


2.9 ±0.6 


7 10~ 3 


1.6 


24 



Ml 


MED 


0.032 ± 


.015 


0.8 


M2 


MED 


0.071 ± 


.027 


0.8 


M3 


MED 


0.013 ± 


.008 


0.5 


M4 


MED 


0.048 ± 


.020 


0.5 


M5 


MED 


0.021 ± 


.016 


0.8 


M6 


MED 


0.044 ± 


.029 


0.8 


Ul 


UHD 


0.15 ± 


.05 


0.8 


U2 


UHD 


0.32 ± 


.09 


0.8 


U3 


UHD 


0.075 ± 


.024 


0.5 


U4 


UHD 


0.23 ± 


.07 


0.5 


U5 


UHD 


0.080 ± 


.032 


0.8 


U6 


UHD 


0.17 ± 


.06 


0.8 



0.8t°j 
2.2 ± 1.0 



.010 
.14 



+ .016 
009 
+ .14 



09 
4? +.49 

+ 1.1 
0.9 



1.2 



5.6 ±2.3 
14± 5 

Q4+0-25 
u -°*-0.18 



3.2 

2.5 
6.5 



+2.2 
1.6 
+ 1.3 
1.2 
+3.0 
-2.7 








0.7 


24 


5 10- 


-3 


0.7 


24 







0.7 


24 


9 10- 


-3 


0.7 


24 







1.6 


24 


3 10- 


-3 


1.6 


24 







0.5 


24 


3 10- 


-2 


0.7 


24 







0.7 


24 


5 10- 


-2 


0.7 


24 







2.4 


35 


3 10" 


-2 


2.4 


35 



a One a errors are shown. The effects of PSF merging (Sect. I3.5[) have not been 
taken into account. 
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Table 3 

Three Dimensional Correlation Lengths 



Fit ID 7 


e 




r c ,o a r c (z cS ) a Model b 








[h- 1 Mpc] 



SI 


1.8 


_1 


2 


1 


07 


9.8±0.7 


9.8±0.7 


U03 


SI 


1.8 


-1 


2 


1 


11 


9.4±0.7 


9.4±0.7 


H05 


SI 


1.8 


-3 





1 


42 


4.3±0.3 


10.4±0.7 


U03 


S2 


1.8 


-1 


2 


1 


07 


12.1±0.9 


12.1±0.9 


U03 


S5 


1.8 


-1 


2 


1 


07 


9.4±0.8 


9.4±0.8 


U03 


SG 


1.8 


-1 


2 


1 


07 


11.8±1.1 


11.8±1.1 


U03 


Ml 


1.8 


-1 


2 





87 




5.8j# 


U03 


Ml 


1.8 


-3 





1 


13 


2 9 +0 - 7 




U03 


M2 


1.8 


-1 


2 





87 






U03 


M5 


1.8 


-1 


2 





87 






U03 


M6 


1.8 


-1 


2 





87 


6-9±l:? 


6.9±1;2 


U03 


Ul 


1.8 


-1 


2 





GO 


H.9lJi 


n-9±l:i 


U03 


Ul 


1.8 


-3 








88 






U03 


U2 


1.8 


-1 


2 





GO 


19±3 


19±3 


U03 


U5 


1.8 


-1 


2 


1 


GO 






U03 


U6 


1.8 


-1 


2 


1 


GO 


12.7^:7 


12.7±a:? 


U03 



i Errors are 1 a. 



b Model designation s- U03 lUeda et~ai] (|2003l ) 
H05 lHasinger et all (|2005l ) 



Table 4 
Estimated cts.agn and Bias 



Fit ID 


ZcS 0"8,AGN 


&AGN 


l0gM ha i o a 








h- 1 M Q 



Without integral constraint. . . 

S5 1.07 1.58±;i| 3.5±0.3 13.6±0.1 

M5 0.87 Q.82±;g 1.7±g;| <13.3 

U5 0.60 lA2tf 2 2 I 6±gj 13.5±0.2 



» 


integral constraint. . 






S6 


1.07 1.63+^ 


3.7±0.3 


13.6±0.1 


M6 


0.87 1.19+;^ 


z -°-i.o 


1 q q+0.3 
iJ -°-0.7 


U6 


0.60 2.0S±;|i 


q O+0.6 


13.9±0.2 



a The error on Afhaio reflects the statistical error 
on &agn only. 
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